Small polarons in dilute gas Bose-Einstein condensates 
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A neutral impurity atom immersed in a dilute Bose-Einstein condensate (BEC) can have a bound 
ground state in which the impurity is self-localized. In this small polaron-like state, the impurity 
distorts the density of the surrounding BEC, thereby creating the self-trapping potential minimum. 
We describe the self-localization in a strong coupling approach. 

PACS numbers: 03.75.Hh,67.40.Yv 



Experimentalists are pursuing the localization and 
transportation of individual atoms in dilute gas Bose- 
Einstein condensates (BECs) 0, . Their motives are 
manifold: The transportation of particles into and out of 
a localized state would realize a quantum-dot-like single 
particle device. The rate at which the localized state re- 
ceives or emits particles could determine the local density 
of states like a scanning tunneling microscope (STM). 
The motion of a localized atom could test superfluid 
dynamics 0, and its acceleration the Unruh effect 0. 
Light resonant with multiple localized particles could it- 
self exhibit localization [5j. The spins of localized par- 
ticles could make up a quantum register of movable qu- 
bits. However, the challenge of localizing a neutral atom 
by means of a steep external potential is daunting 0. 
In this Letter, we propose an alternative strategy: an 
impurity self-localizes within a region smaller than the 
BEC-healing length when the magnitude of the impurity- 
boson scattering length is increased above a critical value. 
Similar to the electron self-localization in a polar crystal 
(forming a small polaron), the BEC-impurity localizes 
because the interaction energy gain (stemming from the 
local distortion of the boson-field) outweighs the kinetic 
energy cost. Observing this phenomenon in cold atoms 
may require a Feshbach resonance (to alter the impurity- 
boson interaction) and impurity creation (either by a Ra- 
man process or by trapping a different atom species), 
but these techniques have been demonstrated p. This 
experiment would create a novel class of self-localized 
particles: polarons with mass comparable to, or possibly 
larger than that of the boson particles. 

In the context of condensed 4 He-fluids, Miller et al. 
remarked on impurity self-localization and mentioned 
the polaron connection. They advocated a perturbation 
treatment (weak-coupling theory) by demonstrating the 
similarity of the perturbed wavefunction with the varia- 
tional one proposed by Feynman and Cohen to include 
'backflow' The weak-coupling description predicts 
that phonon-drag increases the impurity mass from mi 
to m*j = m// [1 — a], where a is a dimensionless coupling 
constant. Its value depends on the BEC density ps, the 
impurity-boson interaction potential V q , the boson mass 
tub, and the excitation energy E^ of the boson modes 



of momentum q: 
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For a dilute BEC and an impurity-boson contact interac- 
tion of scattering length ajs, V q — > \jb = 2irh 2 (l/rriB + 
l/mi)ajB and E^ — > hcsqy/^ + (£<?) 2 - Here, £ is the 
BEC-healing length which depends on the boson-boson 
scattering length cibb, £ = 1 / \f\^p~BOBB ', and cb de- 
notes the BEC velocity of sound, cb — fi/(2mfl(), so that 
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where F(y) = / °° x 5 / [VI + x 2 (xVl + x 2 + yx 2 ) 3 ] . As 
in polaron physics, the effective mass divergence at a = 1 
indicates self-localization, even though the weak-coupling 
description breaks down when a ~ 1 



We describe the self-localized impurity in a strong- 
coupling treatment - similar to the Landau-Pekar de- 
scription of small polarons [lflj - using a product wave- 
function: 

*i,jv(r,xi,x 2 , ...,xjy) ^ x(r)?/>(x 1 )'0(x 2 )...'(/'(x A r), (3) 

where x( r ) represents the impurity wavefunction and 
i> denotes the single particle state occupied by the N 
indistinguishable bosons of position Xj. We substi- 
tute Eq.(PJ into the Hamiltonian eigenvalue equation 
and multiply first by ip* (x2)...ip* (tcn), then by 
X*(r)V ; *( x 2)---V'*( x w)- Integrating the first equation over 
X1X2...XJV, the second over r X2...xjv, and choosing the 
ground state wavefunction to be real valued (e.g. |x| 2 = 
X 2 ) we obtain [llj 
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where Xbb = [47rfi /tob]obb and f/ 3 i s the conden- 
sate field, if = \Nip. If Ein and Bo.jv are the 
ground state energies of N bosons in the presence of 
one or zero impurity atoms, the BEC chemical poten- 
tial is fiB = -E-i.jv — Ei t N—i, and the impurity energy is 
Ei = Ei t N — Eq,n- 

Since the BEC experiences the density of only a single 
impurity, its field may be altered only slightly for suf- 
ficiently weak boson-impurity coupling, ip(r) = \f~fPg + 
S(p(r). Using [ib ~ XbbP%, the corresponding lineariza- 
tion of Eqs. J2J gives 
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where E(, = Ej — XibP% 1S the binding energy. 

As a modified Helmoltz equation, we solve Eq. J5 
terms of the Green function (r) = (An 
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The excess number of BEC-atoms in the impurity region 
is then f dr 2^fpg5<p(r) = — (Xib/Xbb), m agreement 
with [l5j where this number, induced by a general po- 
tential, was determined from thermodynamic arguments. 
The substitution of © into l(5al) results in the wave equa- 
tion of a particle that self-interacts through an attrac- 
tive Yukawa (or screened Coulomb) potential. Exploiting 
the Coulomb analogy, we introduce an effective charge Q 
where Q 2 = [X IyB p B ] a i,B x 2(1 + m B /mi), 
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X(r). 



Incidentally, V mcd (r) = -Q 2 e-l r l/^/|r| is also the BEC- 
mediated interaction experienced by distinguishable par- 
ticles immersed in a BEC, as calculated from perturba- 
tion theory 0, . The Coulomb analogy also suggests 
natural units Eq, the Rydberg energy, and Rq, an effec- 
tive Bohr radius, Q 2 /Rq = h 2 /miR 2 = 2Eq, 
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which set the relevant energy (Eq), time (H/Eq) and 
length (Rq) scales. Note that [ATra 2 B p B }~ 1 is the mean 
free path of an impurity moving among hard-sphere scat- 
terers distributed at the BEC-density. 



The small polaron then corresponds to solutions of 
with negative eigenvalue Ef,. We break transla- 
tional symmetry by hand and solve iteratively for 
an s-wave impurity wavefunction centered on the origin. 
At each iteration step, we solve the eigenvalue problem 
for an impurity particle experiencing an effective poten- 
tial u(r) — J dx |~^ x | e «~~ X 2 ( x ) , in which we sub- 
stitute the impurity density x 2 ( x ) obtained from the 
previous iteration. Spherical symmetry u(r) to a one- 
dimensional integral. Working in natural units iJSJ and 
defining (3 = £/Rq, u(r) reads 
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dr' e- r ' l!i r' X 2 (r') , (9) 



where represents the only density/interaction depen- 
dence that remains, thereby becoming the relevant di- 
mensionless coupling constant, 
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Another candidate, Q 2 /[Hcb], is simply proportional, 
Q 2 /[Hcb] — 2(m,B /mi)f3. In Fig. 1 we show the iter- 
atively obtained x(f )-functions, whereas the inset shows 
the corresponding binding energies (in units of Eq), for 
several /3-values. Thus, this strong-coupling description 
predicts the impurity self-localizes if > 4.7. The deeply 
bound variational wave function (0 > 20) with width 
a = 'i^Jir/2, shown by the bold line of Fig. 1, is remark- 
ably similar to the iterative function [lfij . 

When > 20, the impurity state converges to that of 
a particle bound by a pure Coulomb-self interaction with 
energy E b = -0.316.Bo and extent ^r 2 ) = 4.64i? - 

To understand the interesting 'transition' regime, 
4.7 < < 20, which exhibits an intricate interac- 
tion dependence, we approximate the impurity wave- 
function variationally. The effective impurity equa- 
tion Q is equivalent to minimizing the functional 
E v = T + V/2 Q, where T denotes the kinetic en- 
ergy T = -(h 2 /2mi) J drx(r)V 2 x(r) and V the self- 
interaction energy V = J drx 2 (r)u( r ) , with respect 
to variations of the real-valued normalized wavefunc- 
tion, x( r )- Choosing a Gaussian trial wavefunction, 
x(r) = exp (— |r| 2 /2cr 2 )/(7T(T 2 ) 3 / 4 , the functional, writ- 



ten in natural units, becomes 
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where f(a) — J dr re~ ra e~ r I 2 . Numerical minimiza- 
tion of l|ll|) with respect to a gives a binding energy 
that agrees well with the iterative solution of , shown 
in dotted line in the inset of Fig. 1. The dashed line 
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FIG. 1: Radial wavefunction obtained through the iterative 
procedure for £/R = 4.7,5,10,20,30,40 (from bottom to 
top). In bold black line, the initial Gaussian guess. In the 
inset, the energy of the ground wavefunction vs £ (dots). In 
dotted line, the variational result obtained numerically. The 
expansion for large £/Ro is in dashed line, and the best fit 
E b /E ~ -1/tt + 3R /2£ in solid line. 
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{tub /mi) predicted by the weak (a > 1) 7] and strong 
((3 > 4.7) coupling descriptions. Although neither treat- 
ment should be quantitatively correct, they give compa- 
rable results for 1 < (m B /mi) < 10. 

In an inhomogeneous BEC confined by a trapping po- 
tential Vb (r) , the impurity-free density p° B (r) varies spa- 
tially. Assuming that p° B (r) varies slowly on the scale of 
Ro (i?o|V p B \l Pb « 1)j we describe the self-localized 
impurities (which appear point-like to the BEC) as im- 
mersed in a locally homogeneous superfluid. If the im- 
purities localize on a time scale shorter than the time for 
the impurity to move appreciable (Eq/H » totrap), or 
for the impurities to attract each other (which depends 
on the average impurity density), we can describe the 
subsequent impurity dynamics as that of classical point 
particles subject to an effective potential. This potential 
energy Vj (r) is the sum of the external impurity po- 
tential Vf xt (r), the mean-field energy \jBP B (r), and the 
local binding energy E b \p B (r)\ of Ea. (|12|) . Computing 
Pb{ v ) m t ne Thomas- Fermi approximation, we find 
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where 
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(14) 



Even when Vf xt = and the impurity-BEC interaction is 
repulsive Xjb > - so that the boson mean- field (the first 



FIG. 2: Minimal \/p%aj B /aBB - value for localization of the 
impurity as a function of the mass ratio ms /mi obtained 
from weak coupling (a > 1) Q and strong coupling (Eq. J7J 
and /3 > 4.7) descriptions, dashed and solid line respectively. 



plots the energy obtained by expanding (|11|) for large 
(3, Eb/E ~ — 1/tt + 2//3, showing reasonable agree- 
ment with the iterative solution but slightly overesti- 
mating the minimal /3-value for self-localization. A fit 
that gives better agreement over the whole /3-range is 
E b /E ~ -1/tt + 1.5/(3 (solid line in inset of Fig. 1). 
With Eq = [Xbb /Os]2(to_b / 'raj) (3 2 , E b also equals 



Eb ~ [XBBP B ]^{ m B/mi) 



Z 7T 



(12) 
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reminiscent of the strong coupling energy of traditional 
polarons, proportional to the square of the coupling con- 
stant 0. 

In Fig. 2 we compare the corresponding minimal 
y/pOgajg/aBB-value for localization as a function of which the localization condition, /3 > 4.7, is not fulfilled. 



FIG. 3: Effective potential of the impurity as a function of 
distance to the trap center for (a) attractive and (b) repul- 
sive boson-impurity interaction, using p l g(0)|af B | = 10 -3 , 
\(iib\I clbb = 10 3 and mi = m_g. (c) and (d) have ais > 0, 
but the former has P B (0)\a 3 IB \ = 0.05 and mi = ms, while 
in the latter p a B {Q)\a?i B \ — 10 -3 and m B = 14.5m/ ( 6 Li impu- 
rities in a 87 Kb BEC). The dashed lines show the values at 
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term in the {}-bracket) would expell the impurity from 
the trap center - the binding energy (the other terms 
in the {}-bracket) can give an overall potential that at- 
tracts the impurities to the trap middle. This behavior 
is illustrated in Fig. 3 for typical experimental parame- 
ters for 6 Li impurities in a 87 Rb BEC (m B /mj = 14.5). 
Even if in the true ground state the impurity would hover 
at the edge of the BEC, the self-localization can form 
a metastable state with long tunneling times. In any 
case, Vf xt can keep the impurities within the BEC, and 
Xbi > O-impurities tend to gather in the trap center. 

A question remains regarding the accuracy of the prod- 
uct state © when m B /mj ~ 1, although \E^\ 3> p B 
implies a separation of time scales that justifies the lack 
of impurity-BEC correlations. For \E\>\ < p B , a more 
sophisticated description could be useful: we expect the 
above results to serve as a benchmark for future calcula- 
tions. 

When is the linearization of tp justified? With x(r) = 
e~ r l 2a /(na 2 ) 3 / 4 , the ratio of the peak-value of the 
impurity-induced fluctuation, 5tp(r = 0), to its spatial 
average, equals 5ip(r = 0)/y/fJg = — (l/y / 7r)(l + 

rriB /fni)ifliB / c) / (c /[\/2£]) . Assuming a deeply bound 
polaron, /(c/[\/2£]) ~ 1 an d er ~ 3y^n~/2R 0: the condi- 
tion \5ip(r = 0)/y/^g\ < 1/10 takes the form 
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(15) 



A large increase in aiB above the critical value for self- 
localization could collapse the system when the lineariza- 
tion condition Sp B /p B is violated, as found in ^l|. We 
speculate that in this regime Xjb > impurities could 
'phase separate', creating a hole in the BEC-density. 

In addition, the self-localization condition, (3 > 4.7, 
gives a lower bound to Ipb^bI anc ^ 
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defines the regime in which the linearization assump- 
tion holds and the strong-coupling description predicts 
self-localization. Furthermore, the (ai B /a BB ) ratio is 
subject to conditions stemming from the validity of the 
contact description of the impurity-boson interactions, 
\p B cij B \ << 1, and from the outer-most inequalities of 
(|16fl . giving, respectively 
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These conditions may require a Feshbach resonance, but 
this can be achieved with existing technology. 

Time of flight measurements or diffraction of light reso- 
nant with impurities can detect small polaron formation. 
In the former case, the tightly bound impurity wave func- 
tion can expand faster and further than if the impurity 



were not self-bound ^3- l n the latter case, the opening 
angle 9 of the cone in which light of momentum k is scat- 
tered coherently, 2sin(0/2) ~ 1 / [k yjlr 1 )] , for k ~ R^ 1 
widens abruptly when the impurity self-localizes. 



In summary, we have pointed out that a neutral im- 
purity atom immersed in a homogeneous (or large) BEC 
can self-localize in a region smaller than the BEC-healing 
length. In a strong-coupling description with BEC lin- 
earization, the localizing BEC-distortion gives rise to an 
attractive self-interaction with a spatial dependence iden- 
tical to the BEC mediated impurity-impurity interac- 
tion. Roughly, binding occurs when the range of the 
self-interaction range exceeds the extent of the bound 
impurity - more precisely, when (3 > 4.7, a condition 
that can be fulfilled experimentally. Using a variational 
Gaussian impurity wavefunction, we construct an ana- 
lytical approximation to the binding energy from which 
we obtained the effective potential energy experienced by 
self- localized impurities in a trapped BEC. 
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